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Abstract 

A two-parameter family of complexity measures C^ a '^ based on the Renyi entropies is intro- 
duced and characterized by a detailed study of its mathematical properties. This family is the 
generalization of a continuous version of the LMC complexity, which is recovered for a = 1 and 
(3 = 2. These complexity measures are obtained by multiplying two quantities bringing global 
information on the probability distribution defining the system. When one of the parameters, a 
or (3, goes to infinity, one of the global factors becomes a local factor. For this special case, the 
complexity is calculated on different quantum systems: H-atom, harmonic oscillator and square 
well. 

PACS numbers: 31.15.-p, 05.30.-d, 89.75.Fb. 
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1. INTRODUCTION 



The study of statistical measures in physical systems, and in particular in quantum 
systems, has a role of growing importance. So, information entropies and statistical com- 
plexities have been calculated on different atomic systems l|, |2J. In particular, the so called 
LMC complexity 0, 4] has been computed in the position and momentum spaces for the 
density functions of the Hydrogen-like atoms and the quantum isotropic harmonic oscillator 
. It has been found that the minimum values of that statistical measure is taken on 
the quantum states with the highest orbital angular momentum, just those wave functions 
that correspond to the Bohr-like orbits in the pre-quantum image. 

Many LMC-like statistical complexities are defined as a product of two factors, one of 
them measuring the broadening of the distribution that defines the system and the other one 
quantifying the narrowness of it. Both factors are global magnitudes that can be calculated 
by integrating over the whole support of the distribution. 

Shannon information [?J is an adequate indicator to grasp the spreading of a distribution 
and thus it is employed as a basic ingredient of the first factor of complexity measures. 
Concretely, it plays an important role in the original LMC statistical complexity in which 
the second factor, the so called disequilibrium 3J, is the square distance to the equiproba- 
bility distribution. Different generalizations of the Shannon information that depend on a 
parameter can be found in the literature BQQ. For stance, the Renyi entropy, that eau 
be related with relevant physical magnitudes in atomic physics [HI, H2j , is a good candidate 
to be used in this purpose. 

Hence, in this work, we undertake the generalization of the second factor of the LMC 
complexity by means of the Renyi entropy. We apply the same procedure than in [jj]] where 
the first factor was generalized. This is presented in Section 2. Some properties of this new 
two-parameter-dependent complexity measure are indicated in Section 3. Strikingly, when 
one of the parameters tends to infinity, the asymptotic limit of this measure becomes the 
product of a global quantity by a local one. In Section 4, the analysis and calculation of 
this special global/local product case of the generalized complexity measure are performed 
for the H-atom, the quantum harmonic oscillator and the square well. Last section includes 
the conclusions. 
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2. GENERALIZED STATISTICAL COMPLEXITY MEASURE cf ,f)) 

Let us consider a D-dimensional density function /(r), (with /(r) nonnegative and 
J" f(r)dr = 1). The Renyi entropy of order a of the density function / is given by 



R (a) 



—^—Inflf^rdr, for < a < oo, 1, (1) 
1 - Oi J 



the limit a — > 1 we have Rf — > 5 / = — J /(r) ln/(r)dr, with Sj the Shannon entropy of 
/, and in the limit a — > oo then ir^ — > — In ||/||oo, where H/Hoo = sup r f(r) represents the 



where r stands for n, ...,ro. From the above definition, it is straightforward to see that in 

maximum reached by / over its whole support. 

The importance of Renyi entropies comes from the fact that, for atoms and molecules, 
density functionals (kinetic energy, exchange energy and classical Coulomb repulsion energy) 
can be expanded in terms of the local homogeneous functionals exp (1 — a)R^ {ill. \\\. In 
particular, it is well known that for a = 5/3, a = 4/3 and a = 2, the local density 
approximations are related to the kinetic and exchange energies and the average of the 
density, respectively. 

A continuous version Q] of the measure of complexity Cf, the so-called LMC complexity 
introduced in [3[, is defined by 

Cf = H f Q f , with H f = e s f and Q f = e~ R ? . (2) 

When the Shannon entropy of the statistical complexity Cf is replaced with the Renyi 
entropy of order a, we obtain the generalized statistical measure of complexity, Cj , which 
has been defined by {3| 

Cf ] = Hf ] Q f , with Hf ] = e R[ f a \ (3) 



and tends to Cf in the limit a — > 1. 

Now we can substitute in a symmetric way the i?j ingredient of the above complexities 
for the Renyi entropy of order /3, which allows us to obtain a (a, (3)— dependent measure of 

7 



complexity, Cl Q,/3 \ which is defined by 



C^^e< } - R f\ 0<a,/3<oo. (4) 

bo, we recover Cf 2) = Cf and Cf a) = Cf of Refs. , respectively. This type of 

generalization based on Renyi entropies differences was suggested in [l^] after the work of 



Varga and Pipek 



151 ] in this same line of thought. 



3. PROPERTIES OF cf' P) 



Now, having as guideline Ref. 



4|, we proceed to present some mathematical properties 



of this new generalized statistical complexity measure. 

A. Inversion symmetry 

It is straightforward to check that 



7 w / 

and then Ci a ' a) = 1. 



B. Monotonicity and universal bound 

Taking into account that the Renyi entropy is a nonincreasing function of a, it can be 
easily proved that 

(i) Cf' P) > 1 if a < (3 and Cf' P) < 1 if a > 0. 

(ii) Ci a is a nonincreasing function of a for a fixed (3 and an increasing function of j3 
for a fixed a. 

(iii) The lower (upper) bound 1 is reached for a < (3 (a > (3) for the rectangular density 
function. It is a universal bound independent of a and /3, as it is shown in Sec. 13 Fl 



C. Invariance under translations and rescaling transformations 

Cj"'^ is invariant under scaling transformations, i. e. for fx = A D /(Ar), then C^'^ = 
(jy>p) _ Also, it is invariant under translations. Hence, in general, let us consider a scaling 
transformation and a translation parameterized by (a,b), respectively, of the form 

/ab(r) = a D /[a(r-b)], (6) 

where the distribution function / is normalized, j f(r)dr= 1. 
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The Renyi entropy of order a is transforming as 

= In J (a D f[a(r - b)]) a dr = R {a) - Din a, (7) 

where a change of variable y = a(r — b) was applied. Note that there is no dependence on 
the parameter b. Therefore the new complexity measure (jlj) is invariant under this 

transformation. 

D. Invariance under replication 

Take n copies / m (r),?7i — 1, ...,n of the distribution function /(r), 

/ m (r) = n^ 2 -7K /2 (r - b m )], 1 < m < n, (8) 

where the support of each f m (r) is centered at the point h m and the supports are disjoint. 
We can immediately obtain that J f m (r)dr = 1/n. Then we can define the distribution 



function 



q(r) = J2Ur), (9) 

m=l 

that is normalized, J q(r)dr = 1. From Eq. (jSJ) we can easily calculate that 

J fZ(r)dr = n^ D / 2 ~ a J T(y)rf y . (10) 
As the replicas are supported on disjoint sets, we have 

[ r m {v)dv = n^ D i 2 ^ I r(y)dy. (ii) 

m=l ^ 

Then the transformation of the Renyi entropy of order a is 

Rf) = R^ - (D/2- l)lnn. (12) 
Hence the complexity measure Cq a '^ (jl]) is replica invariant. 

E. Near-continuity 

Take two distribution functions /(r) and g(r) defined on the set M in the D dimensional 
space considered. Let 5 be a positive real number. The functions /(r) and g(r) are 5- 
neighboring functions on M, if the Lebesgue measure of the points r e M satisfying |/(r) — 



g(r)\ > 5 is zero. A functional T of the distribution functions is near-continuous if for any 
e > there exists 5(e) > such that for any 5-neighboring functions /(r) and g(r) on T 
then \T(f)-T(g)\<s. 
Take the function 

— if |r| < 1 

•| < B , (13) 



CD 

^31) i/ K M < B 
otherwise 



where 



2tt d / 2 



c ^ = Drw^)' (14) 

where B > 1 and 1 > 5 > 0. As cd is the volume of a unity .D-dimensional sphere then g is 
normalized to 1. One can easily calculate the Renyi entropy of order a, 



R (a) 



1 — a 



In 



:i-sr+ 



5° 



(b d - 1; 



a-l 



+ In cd . 



(15) 



From here, we obtain C, 
density function 



(«./») 
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< a,j3 < oo. Consider now the rectangular 



The Renyi entropy of order a 



— if |r| < 1 

CD 

otherwise 



=\nc D 



(16) 



(17) 



does not depend on the a. Consequently, we are led to the result 



n(a,0) _ i 
X 



(18) 



Note that g$ b and \ are 5-neighboring functions for < 5 < 5 < 1 and 



Ums^oRW = R^ = In c D . 



(19) 



Therefore in the limit 5^0 the complexity measure takes the same value: 



lim s ^ & a '® = Ci a '® = 1. 

y a, 



(20) 



The importance of a bounded support to obtain this result 
to be done in a future work, such as it was suggested in 



; deserves some longer explanation 
1 
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F. The extremal complexity 



The extremal complexity is reached for the rectangular function (Eqs. (fT6|) and (ITS]) ). It 
can be proved following J|. The function / is taken as a sum of rectangular functions XE k 
defined on disjoint sets E k , k = 1, ...,n with Lebesgue measure fi^ 

n 

f(r) = Y,hXE k . (21) 
fc=i 

Its integrals can be easily calculated: 

/n 
f a dr = K»k ■ (22) 
k=i 



The logarithm of complexity measure C^ a ^ has the form 



W . (23) 



We seek the extremum of the logarithm of complexity measure C( a '@> under the normalization 
condition 

n 

J2 X ^k = l- (24) 
k=i 

The variation with respect to A& and fik after straighforward manipulation leads the equa- 
tions 



1=1 X l 



for all k — 1, n. As has the same value for all k, then / is a rectangular function. 
Moreover, with the help of the near-continuity property as explained in 4|, it can also 
be argued that the rectangular distribution is the only distribution reaching the extremal 
complexity. 



G. The cases C'j Q ' oo) and cj a ' 0) 

As before explained, the new complexity measure has a completely different behavior for 
a < (3 and a > /3. In the first case, there is a lower bound, and in the second case there is an 
upper bound. Both are universal, i. e. the bound is equal to 1 for any (finite and non zero) 
value of a or /3, and this universal bound is reached at the rectangular density function. 
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This is a consequence of the fact that the Renyi entropy is independent of the parameter a 
or (3, then the new complexity measure has the value of 1 for any a or (3. 

Finally, let us remark that when j3 goes to infinity a special case of the complexity measure 
in terms of a local quantity (||/||oo = su Prf( r )) is obtained: C , j a '°°' ) = e R f ||/||oo- To prove 



this property it is sufficient to take into account that linip^oo ( f /(r) p c/r) 1 ^ p = sup r f(r) 16]. 
This complexity measure verifies Cf )0o) > 1 for all a (and C [ p a) < 1). We can also mention 



j -L LKJL an lx ^aiiu V_y f 



that cl a ' ^ — > for finite a and /3 — > 



4. CALCULATION OF cf ,oo) ON DIFFERENT QUANTUM SYSTEMS 

Among the different statistical indicators that have been defined as a product of two en- 
tropic terms, the generalized complexity C^ a '^ is also obtained by multiplying two factors, 
each one of them bringing global information on the distribution /. Recall, however, that the 
limit C , j a '°°' ) combines global information on the distribution /, just the part corresponding 
to e f , and local information coming from only a specific point of the space where the dis- 
tribution / is supported, indeed the maximum of the density, ||/||oo- Equivalently, the limit 
version C , j°°' /3 ' > presents a symmetrical behavior respect to ci 01 ' 00 ^ by considering the prop- 
erty C^'^Cj 7 ' 00 ^ = 1. Here we undertake the calculation of C , j a '°°' ) for three paradigmatic 
quantum systems: the if -atom, the harmonic oscillator and the square well. 



A. .ff-Atom 



The atomic state of the if -atom determined by the quantum numbers (n, l,m) in position 
space (r = (r, Q), with r the radial distance and Q the solid angle) is given by the non- 
relativistic wave function 

*n,i,m(r) = Rn,i( r ) Y i,M , (26) 

where R n ,i{r) is the radial part and Y^ m (Q) is the spherical harmonic. The radial part is 
expressed as [171 ] 



ii- 



{n-l-l)\ 
(n + l)\ 



1/2 / 2r W 



n 



-- T 2l+l 
e " L n-l-l 



2r 



ii 



(27) 



L@(t) being the associated Laguerre polynomials. Atomic units are used in this section. 
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The same hydrogenic atomic state in momentum space (p = (p, Cl) , with the momentum 
modulus p and the solid angle Cl) is given by the wave function 

Khmip) = RnM Y lttn (fl) , (28) 

where the radial part R n ,i{p) is expressed now as [l 
'2 (n-l-l)\ 



Rn,l(P) 



1/2 ; ; / 2 2 i ' 

(nV+l) ,+2 "^InV + l 1 ' 



_7T (n + Z)! 

with the Gegenbauer polynomials G^(t). 

Taking the former expressions, the probability density in position and momentum spaces, 

p(r) = | ^, m (r) | 2 , 7(P) = I *n,z, m (p) | 2 , (30) 

can be explicitly calculated. From these densities, we compute Cf* , taking into account 
that for these cases / = p(f) or / = j(p), respectively. 

In Figure 1, the value of the generalized complexity in position space, Cr a '°°\ is shown 
for n — 15 and I = 5, 10, 14 versus \m\ with a = 0.5 (Fig. 1(a)) and a = 2 (Fig. 1(b)). The 
same calculation in momentum space, Cp ' \ can be seen in Figs. 2(a) and 2(b), for the 
cases a = 0.5 and a = 2, respectively. Note that the minimum of the generalized complexity 
corresponds just to the highest I for a given n in both position and momentum spaces. This 
property is independent of the parameter a. 



B. Harmonic oscillator 



Let us consider a particle under the action of the potential energy V(r) = A 2 r 2 /2, where 
A is a positive real constant expressing the potential strength. The three-dimensional non- 
relativistic wave functions of this system in position space (r = (r,d)) are: 



(31) 



where Rn,i(r) is the radial part and Yi jm (fl) is the spherical harmonic of the quantum state 
determined by the quantum numbers (n,l,m). The radial part is expressed as [17 ] 

1/2 



RnAr) 



2 n\ A'+ 3 / 2 
_r(n + Z + 3/2) 



re 2 



(32) 
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where L^(t) are the associated Laguerre polynomials. The levels of energy are given by 

E n .i = X(2n + I + 3/2) = A(e n ,z + 3/2), (33) 

where n = 0, 1, 2, • • • and / = 0, 1, 2, • • • . Let us observe that e n j = 2n + I. Thus, different 
pairs of (n, I) can give the same e n j, and then the same energy E n ^. 
The wave functions in momentum space (p = (p, Cl)) are: 

Khm(p) = RnM Yi, m (Ci) , (34) 
where the radial part R n j(p) is now given by the expression 



Rn,l(P) 



V 



1 e-fe LL +1/2 (P 2 /A) • (35) 



_r> + Z + 3/2)_ 

Taking the former expressions, the probability density in position and momentum spaces, 

P\{r) = | ^n,l,m(r) T , 7A(P) = I ^nJ,m(p) ? , (36) 

can be explicitly calculated. From these densities, the generalized statistical complexity is 
computed. The subindex A can be dropped because this indicator, C^ a '^\ is independent 
of the potential strength, A, due to its invariance under scaling transformation. As a con- 
sequence of this property, it is also found that this magnitude is the same in both position 
and momentum spaces, Ci a '^ = Cp 01 '® . 

In Fig. 3, Cr 01 (or Cp a '°°^) is plotted as function of the modulus of the third component 
m, —I < m < I, of the orbital angular momentum I for different I values with a fixed energy, 
e n j = 15, when a = 0.5 in Fig. 3(a) and a = 2 in Fig. 3(b). It can be observed that 
Cr a '°°^ splits again in different sets of discrete points. Note that the values associated with 
the highest I (I = 15) give the minimum values of Ci a '°°\ 



C. Square well 



The eigenstates of the energy in the quantum infinite square well in position space for 
a particle in a box, that is confined in the one-dimensional interval [0,L], are given by the 
wave functions [191 ] 



<Pk{%) 



kirx 



1,2, 



(37) 
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The probability density of the kth excited state is 



p(x) =| ip k (x) 



(38) 



that gives a maximum of 2 when L is considered as the natural length unit in this problem. 
The other factor necessary to obtain the generalized statistical complexity, Cr a '°°\ gives 



Then, Cr = 2g(a) with g(a) — e f . We conclude that the statistical complexity 
is degenerated for all the energy eigenstates of the quantum infinite square well. Its value 
can be computed as a function of a. It takes 2 for a = and decays monotonically to 1 
when a — > oo. In the general case of a particle in a d-dimensional box of width L in each 
dimension, it can be also verified that complexity is degenerated for all its energy eigenstates 
with a constant value given by Cr a = (2 g(a)) d . 

5. SUMMARY 

The generalized complexity measure defined here provides a family of complexity mea- 
sures. We have performed a detailed mathematical characterization of its properties. As 
usual, these complexities have been defined by multiplying two factors, each one bringing 
global information on the probability distribution defining the system. The whole family 
is identified by two parameters, a and (3. For the special case of j3 going to infinity, it is 
remarkable that the new complexity measure is the product of a global quantity by a local 
information of the density distribution. Then we have carried out the calculation of C'j"' 00 ' 1 
for different quantum systems: H-atom, harmonic oscillator and square well. We have found 
that the behavior of the complexity for these systems in this specific case of global/local 
terms product is similar to that displayed by it in the general case of global/global terms 
product. 




(39) 
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FIG. 1: Generalized statistical complexity in position ' , vs. I to I for different I values 

when n = 15 in the hydrogen atom, (a) a = 0.5 and (b) a = 2. All values are in atomic units. 
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FIG. 2: Generalized statistical complexity in momentum space, Cp a,0 °\ vs. \m\ for different Z 



values when n = 15 in the hydrogen atom, (a) a = 0.5 and a = 2. All values are in atomic units. 
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FIG. 3: Generalized statistical complexity in position space, Cj a '°°\ and momentum space, Cp"' 00 "*, 
vs. \m\ for the energy c n ^i — 15 in the quantum isotropic hsirmonic oscillator for (a) o. — 0.5 cincl 
(b) a = 2. Recall that c£* ,oo) = C$ vx>) . All values are in atomic units. 
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